Ultrafast surface carrier dynamics in the topological insulator Bi 2 Te 3 
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We discuss the ultrafast evolution of the surface electronic structure of the topological insulator 
Bi2Te3 following a femtosecond laser excitation. Using time and angle resolved photoelectron spec- 
troscopy, we provide a direct real-time visualisation of the transient carrier population of both the 
surface states and the bulk conduction band. We find that the thermalization of the surface states is 
initially determined by interband scattering from the bulk conduction band, lasting for about 0.5 ps; 
subsequently, few ps are necessary for the Dirac cone non-equilibrium electrons to recover a Fermi- 
Dirac distribution, while their relaxation extends over more than 10 ps. The surface sensitivity of 
our measurements makes it possible to estimate the range of the bulk-surface interband scattering 
channel, indicating that the process is effective over a distance of 5 nm or less. This establishes a 
correlation between the nanoscale thickness of the bulk charge reservoir and the evolution of the 
ultrafast carrier dynamics in the surface Dirac cone. 



PACS numbers: 78.47.J-; 79.60.-i; 73.20.-r. 

Topological insulators are new states of quantum mat- 
ter that are attracting a lot of interest due to their unique 
transport properties [TJ [5] . A distinctive feature of these 
materials is the fact that they have conducting states 
at their edge or surface which are protected by time- 
reversal symmetry: as a consequence, they cannot be 
adiabatically connected to the bulk which is instead in- 
sulating [3l 0] . Examples of three-dimensional topologi- 
cal insulators that present these properties - with their 
metallic surface bands constituting a Dirac cone - are 
Bi 2 Se3, Sb2Te3 and Bi 2 Te3 [5]. These properties make 
them natural candidates for a new generation of elec- 
tronic devices [6], but for such a development to be real- 
istic a description of the dynamical response of these ma- 
terials and a detailed knowledge of their conduction elec- 
tronic states are needed. Furthermore, since the Dirac 
cone is a surface feature, it is essential to determine its 
interband scattering range, i.e. the distance over which 
the surface bands effectively exchange carriers with the 
bulk. 

In this letter, we show how pump-probe ultrafast An- 
gle Resolved PhotoElectron Spectroscopy (ARPES) can 
be used to answer these questions. ARPES is playing a 
key role in the study of topological insulators, thanks to 
its ability to map electronic states in k-space and to its 
surface sensitivity [7HTP]. The same unique advantages, 
coupled with the femtosecond resolution achievable with 
laser sources in pump-probe configuration, make time- 
resolved ARPES the technique of choice to explore the 
transient electronic structure and the evolution of out-of- 
equilibrium topological insulators. In the case of Bi2Te3, 
we found that the bulk conduction band acts as a charge 



reservoir that reaches a quasi-equilibrium condition with 
the surface bands in about 0.5 ps and allows us to esti- 
mate the ballistic range for electrons travelling within the 
surface in the edge bands (180 nm). The hot electrons in 
the Dirac cone evolve towards a Fermi-Dirac distribution 
over few ps, while their subsequent relaxation process, 
lasting for more than 10 ps, is considerably longer than 
the main bulk recombination processes. Furthermore, 
the comparison of our time resolved ARPES results with 
ultrafast optical reflectivity measurements on the same 
material indicates that in Bi2Te3 the interband scatter- 
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FIG. 1: (a) ARPES image acquired along the TK direction 
in s-polarization with the 6.3 eV laser source, (b) pictorial 
view of the experiment: the infrared laser pump excites elec- 
trons from the sample valence band (orange) into the empty 
conduction band (light blue). The subsequent carrier flow 
produces a transient variation in the charge population in 
the surface Dirac cone, both for empty (blue) and filled (red) 
electronic states, and is probed by the ultraviolet probe. 
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FIG. 2: Time-resolved ARPES sequence obtained after photoexcitation from the pump pulses: the color intensity is presented 
in logarithmic scale to make the signal from the transient electronic states more evident. In the upper left corner, a schematic 
view of the bands involved in the process: the projections of the bulk bands B* and B%, and the surface Dirac cone. In the 
first image (t=-200 fs) the integration regions S, ST, B* and are defined 



ing of carriers between bulk and surface states is effective 
over a distance perpendicular to the surface of 5 nm. This 
can be a relevant piece of information in possible future 
device applications based on thin films of this material, 
since it means that the film thickness can affect the ultra- 
fast response of the Dirac cone, representing a potential 
tuning parameter to control carrier dynamics to and from 
the surface topological bands. 

Single crystals of Bi2Te3 (n-type) were carefully ori- 
ented and cleaved under ultrahigh vacuum conditions in 
the recently commissioned FemtoARPES experimental 
setup [TTJ [T2], where pump-probe photoemission mea- 
surements can be performed with a time resolution of 65 
fs and an energy resolution of 60 meV. Photoelectrons 



were excited with 6.32 eV photons, obtained as the 4th 
harmonic of the 1.58 eV laser which is used to pump the 
transient states in pump-probe configuration. The use of 
a 6.32 eV probe gives access to a limited portion of k- 
space, but sufficient in this case to observe the Dirac cone 
produced by the surface bands. In Fig. [T|a) we present 
a conventional (i.e. without the pump laser) ARPES im- 
age along the FK direction. The results are consistent 
with previous ARPES studies at higher photon energy on 
Bi2Te3 [8]: the surface Dirac cone is clearly visible. The 
probing depth of low energy ARPES can be estimated 
to about 2-3 nm [13 . This optimizes the simultane- 
ous detection of the surface dynamics (in particular the 
Dirac cone) and of the bulk dynamics which is probed 
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by ultrafast optical reflectivity measurements on Bi2Te3 
[14l[T5]. The sibling compound Bi 2 Se3 was also studied 
with ultrafast optical spectroscopy [HI [17], and selec- 
tive sensitivity to the surface was achieved using second 
harmonic reflections [T7j. With fs pump-probe ARPES, 
one has direct access to the transient electron popula- 
tion and its density of states: in this work we focused on 
these aspects, paying particular attention to the sub-ps 
thermalization regime. 

In this study we present results obtained on the VK 
direction, where no hybridization is present between the 
surface Dirac cone and the bulk conduction band states in 
the proximity of the Fermi energy Ep [8] [9]: this made 
it possible to unambiguously discriminate the photoemis- 
sion yield between the empty surface (blue in Fig. [TJb)) 
and bulk (light blue in Fig.[ljb)) bands over an extended 
wavevector-energy window. We explored also other direc- 
tions in fc-space (namely TM) , and no relevant difference 
was found in the observed carrier dynamics, indicating 
that the hexagonal deformation of the Dirac cone present 
in Bi2Te3 does not play a critical role in the ultrafast re- 
sponse of the system. All the data presented here were 
taken at room temperature. 

In Fig. [2] we present the evolution of the ARPES yield 
at different pump-probe time delays. The schematic 
band structure of Bi2Te3 shown in the upper left cor- 
ner of Fig. [2] is consistent with recent theoretical calcu- 
lations [3 [T5], and can be used as a guide to the eye 
to interpret our data. For sake of simplicity, during our 
study we followed the evolution of the electron popula- 
tion of the regions in /c-space shown for the t=-200 fs time 
delay in Fig. [2j one for each projection of the bulk con- 
duction bands B{ and J3 2 , one for the occupied surface 
band S and one for the unoccupied surface band S*. In 
particular, the separation between the conduction band 
B\ and the surface band S* is clear from the experimen- 
tal data (for instance at t=700 fs), so that our choice of 
the integration regions makes it possible to unambigu- 
ously follow the evolution of the different bands. This 
approach allowed us to capture the essential features of 
the ultrafast carrier dynamics in the nanometer region 
close to the surface. 

Following the pump excitation (50 fs, 1.58 eV, 0.15 
mJ/cm 2 ), a transient excited electron population is ob- 
served in the conduction band B\ around 0.6 eV above 
Ep (t=100 fs). This transient population relaxes quickly 
(in about 300 fs), triggering a cascade of recombination 
processes, and acting as a source of electrons for the bulk 
conduction band B 2 and for the surface bands. This is 
well visible also in the movie provided with the supple- 
mentary material. The images in Fig. [2] provide a di- 
rect visualisation of the electron dynamics, but in order 
to get more quantitive information, we also provide in 
Fig. [3] the time evolution of the populations of S, S* , B^ 
and £?2 , as well as of the overall excited bulk population 
B*(t) = Bt(t)+B%(t). 



The ultrafast carrier dynamics between the surface (S, 
S*) and the bulk (B*) can be described using a stan- 
dard rate equation model for semiconductors [HI [20] 
to fit the population curves S*(t), S(t) and to identify 
the main charge flow processes to and from the sur- 
face bands: for instance, femtosecond carrier dynam- 
ics could be successfully studied at the surface of III-V 
semiconductors, as clearly demonstrated for the surface 
resonance of InP(100) [3T]. Nevertheless, in the case of 
Bi 2 Te3 one is not confronted with weakly dispersing sur- 
face resonances, normally found in standard semiconduc- 
tors: on the contrary, like for all topological insulators the 
surface bands are metallic, strongly dispersing and cross 
the whole bulk energy gap. This makes it necessary to 
adopt an alternative approach to study their relaxation, 
rather than introducing in the rate equations recombi- 
nation terms proportional to the product S(t)S*(t)) as 
normally done for the description of semiconductor sur- 
face resonances. As we will see later (Fig. [4J , the detailed 
study of the evolution of the photoelectron Energy Dis- 
tribution Curves (EDC) of the Dirac cone can represent 
such an alternative approach. 

Based on the assumption that B\ acts as a source of 
excess electrons that adds up to direct excitation for B 2 
and S* with relative weights a and j3, respectively, the 
main carrier interband scattering processes can be de- 
scribed by the following equations: 
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the flow of charge to S* (t) is described by direct ex- 
citation from the laser pulse Gs*, by a flow of electrons 
from the B\ (t) projection of the bulk conduction band 
and by the recombination processes for excess electrons 
described by T m . Similarly, the flow of charge to S(t) is 
described by direct excitation from the laser pulse Gs 
(gaussian with the same FWHM), by a flow of holes 
from the bulk valence band (B to t — B(t)) with charac- 
teristic time Th, and by the recombination processes for 
excess holes (S to t — S(t j) described by tdi- can 
be also described by direct excitation, by interband scat- 
tering from Bl, and by a decay time r 2 . We used the 
solutions of these coupled differential equations as fit- 
ting functions for our population curves Bl(t), B%(t), 
S*(t) and S(t) , obtaining as a result the curves shown as 
solid lines in Fig. [3] for the following values of the fitting 
parameters: FWHM of 76 fs for G u G 2 , G*s*, and G s ; 
a = 0.18; 13 = 0.12; n = O.Mps; r 2 = l.Sps; t d1 = 2.1ps; 
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r h = Q.lips; t D2 = 5.7ps. 
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FIG. 3: Ultrafast evolution of the populations of S, S* , B{ 
and £>2- The fitting curves (solid lines) were obtained as 
explained in the text. 



There are obviously many different phenomena that 
contribute to the different characteristic times presented 
here (electron-phonon scattering, carrier trapping at de- 
fects, etc.), both for surface and for bulk bands. We 
adopted a simplified description (surface vs bulk) because 
our main interest here is to determine the scattering time 
between the Dirac cone and the bulk bands: as we have 
seen, the fast decay of B\ is strictly related to the time 
needed for both S and S* to reach the maximum pop- 
ulation change (Fig. [3]). S(t) and S*(t) present never- 
theless a slightly asymmetric behavior while reaching a 
thermalization regime, and also during the relaxation of 
the excess carriers, as indicated by the different values 
of Tux and T£>2- Electron- hole recombination between S 
and S* gives of course a symmetric decay term, conse- 
quently the asymmetry must be due to the presence of 
decay terms which are markedly different between elec- 
trons and holes (like scattering to defects). The asym- 
metric flow of electrons and holes on the ps and longer 
time scale from surface states is frequently found in many 
other semiconducting materials, and can be for instance 
simply related to the inequivalent position of the sur- 
face bands with respect to the valence and conduction 
bands [22] : we shall not discuss this here, but we found 
that longer scans show that after about 15 ps both S* (t) 
and S(t) are back to their values before excitation. Fur- 
thermore, the very similar values of tjji and Ti indicate 
that S*(t) and B 2 (t) present a parallel evolution: even 
though we didn't explicitely insert a scattering term be- 
tween them in the rate equation, this similarity indicates 
that an effective carrier exchange takes place between 
the two bands. After completion of our work, we became 
aware of a time-resolved ARPES experiment on p-type 
Bi2Se3 |23) . where the authors also found a transient sur- 
face state population persisting for more than 10 ps, due 



to the reservoir role played by the bottom of the conduc- 
tion band. The fact that different materials with different 
doping levels show this effect suggests that a slow relax- 
ation of surface hot electrons may be a general property 
of topological insulators. 

Following the temporal evolution of the populations 
gives a good description of the charge flow, but more pre- 
cise information can be obtained on the electronic struc- 
ture of the photoexcited surface states by looking at the 
details of the photoelectron EDC's. In Fig. [3] we present 
the EDC's extracted from the surface bands (the momen- 
tum integration regions are visible in the inset) at some 
representative time delays: these curves can give a direct 
description of the Dirac cone during its thermalization 
process. One can see that for the longer delays, namely 
2, 3 and -0.2 ps after the pump-pulse, the EDC's can be 
well fitted with a Fermi-Dirac distribution (the spectra at 
2 and 3 ps give an electronic temperature of 792 K and 
619 K, respectively). For shorter time delays, i.e. 0.1, 
0.2 and 0.4 ps, the EDC's cannot be fitted with a Fermi- 
Dirac distribution, due to a strong presence of electrons 
at high energy, indicating that the system is not ther- 
malized yet: at these delays the population of S* is still 
increasing due to the flow of charge from the bulk bands. 
This confirms that the dynamics of the surface Dirac cone 
is somehow delayed with respect to the bulk conduction 
band dynamics. Furthermore, if one multiplies the ther- 
malization time (0.5 ps) by the Fermi velocity estimated 
from the ARPES data (0.36 m//xs), one finds that the 
scattering distance for electrons travelling parallel to the 
surface is of the order 180 nm: this sets a length scale for 
estimating the ballistic range for transport in the topo- 
logical bands. 
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FIG. 4: Energy Distribution Curves along the surface state 
S+S* at different pump-probe delays. For the spectra at -0.2, 
2 and 3 ps we also present the corresponding fits with a Fermi- 
Dirac distribution. The photoemission yield is presented in 
logarithmic scale. In the inset, the corresponding momentum 
integration area is shown. 
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We should now discuss some implications of these re- 
sults, also comparing our pump-probe ARPES data with 
ultrafast reflectivity in the optical range on Bi2Te3. The 
two techniques are very complementary in terms of sur- 
face sensitivity: while the probing depth of low energy 
ARPES is about 2-3 nm [T3], optical reflectivity is sen- 
sitive to more than 10 nm [T3]. One should also keep in 
mind that while our B* is a signal given only by the con- 
duction band population, optical reflectivity probes the 
joint density of states from both the valence and the con- 
duction band. Ultrafast reflectivity measurements per- 
formed by ourselves on the same single crystals studied 
with photoemission as well as by other authors on thin 
films of Bi2Te3 14J at comparable fiuences make it pos- 
sible to evaluate the carrier density at different depth 
values z from the surface (see Fig. 4 in [2]). A fast de- 
cay time can be identified in the carrier density evaluated 
from reflectivity data [14], which has been attributed to 
electron-phonon scattering (0.5 ps) and to carrier diffu- 
sion (more than 1 ps). It should be noted that even 
from reflectivity, the sub-ps carrier decay is estimated to 
be more pronounced at z=0 than at z=10 nm (Fig. 4 
in [14|). The comparison with our data shows that in 
B* (t) the sub-ps decay is much more pronounced due to 
the charge flow to S* , with r e = 0.35 ps - mainly because 
ARPES is more surface sensitive. This makes it possible 
to obtain a quantitative estimate of the distance from the 
surface over which scattering for excess electrons from the 
conduction band to the surface bands is effective. If one 
assumes an exponential attenuation of the photoemission 
yield I(z), I{z) = Ioexp(—z/fx), one finds that for a sur- 
face signal to be comparable to a bulk one, it must be 
originated from a slab of thickness Az comparable to fi 
(2 nm); furthermore, it must certainly be less than the 
sensitivity of optical reflectivity (10 nm). We can thus 
conclude that in Bi 2 Te3 the scattering of hot electrons 
from the conduction band to the Dirac cone is confined 
to a surface layer with a thickness of less than 5 nm. 
This means that altering the thickness of the subsurface 
region on the nanometer scale will affect the transient 
carrier dynamics in the Dirac cone, providing an impor- 
tant element and a potential tailoring parameter for pos- 
sible technological applications of Bi2Te3 and of other 
topological insulators. In particular, the knowledge of 
the carrier scattering between the surface bands and the 
bulk is critically important to extend to topological insu- 
lators some technological possibilities demonstrated for 
truly bidimensional systems like graphene [23], related 
to ultrafast population inversion in the Dirac cone. 

In conclusion, we studied the transient electronic dy- 
namics at the surface and subsurface region of the topo- 
logical insulator Bi 2 Te3. Using time resolved ARPES, we 
provided a direct visualisation of the excess carrier pop- 



ulation and of its evolution. We found that the ultrafast 
dynamics of the carriers in the surface Dirac cone is de- 
layed with respect to the bulk, because the bulk acts as a 
reservoir that keeps providing a relevant charge flow for 
more than 0.5 ps, and because the subsequent relaxation 
phase (more than 10 ps), governed by charge diffusion 
and electron-phonon coupling, is less efficient than the 
typical bulk recombination. The surface sensitivity of our 
measurements make it possible to provide a quantitative 
estimate of the distance over which the bulk-surface in- 
terband carrier scattering is effective (about 5 nm). The 
results presented here for Bi 2 Te 3 appear of general in- 
terest also for other three-dimensional topological insu- 
lators, in particular for possible nanoscale technological 
applications of these materials. 
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